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ABSTRACT 

Analyses of very accurate CoRoT space photometry, past Johnson V photoelectric photo- 
metry and high-resolution echelle spectra led to the determination of improved and consistent 
fundamental stellar properties of both components of AU Mon. We derived new, accurate 
ephemerides for both the orbital motion (with a period of 11"? 11 3) and the long-term, over- 
all brightness variation (with a period of 416'?9) of this strongly interacting Be -i- G semi- 
detached binary. It is shown that this long-term variation must be due to attenuation of the 
total light by some variable circumbinary material. We derived the binary mass ratio Mg/Mq 
= 0. 17+0.03 based on the assumption that the G-type secondary fills its Roche lobe and rotates 
synchronously. Using this value of the mass ratio as well as the radial velocities of the G-star, 
we obtained a consistent light curve model and improved estimates of the stellar masses, radii, 
luminosities and effective temperatures. We demonstrate that the observed lines of the B-type 
primary may not be of photospheric origin. We also discover rapid and periodic light changes 
visible in the high-quality residual CoRoT light curves. AU Mon is put into perspective by a 
comparison with known binaries exhibiting long-term cyclic light changes. 

Key words: stars: binaries: general - stars: binaries: eclipsing - stars: emission-line. Be - 
accretion - stars: individual: AU Mon 



* Based on photometry collected by the CoRoT space mission as well as 
spectroscopy obtained with the FEROS spectrograph attached to the 2.2- 
m telescope at ESO, La Silla, Chile, under the ESO Large Programme 
LP178.D-0361, and with the SOPHIE spectrograph of the Observatoire de 
Haute-Provence (France). The CoRoT space mission was developed and is 
operated by the French space agency CNES, with participation of ESA's 
RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and 



1 INTRODUCTION 

AU Monocerotis was selected as one of the few known bi- 
nary targets in the asteroseismology field of the French-European 

Spain. Based on observations collected at the Centro Astronomico Hispano 
Aleman (CAHA) at Calar Alto, operated jointly by the Max-Planck Institut 
fiir Astronomic and the Instituto de Astrofi'sica de Andalucfa (CSIC). 
t E-mail: maarten.desmet@ster.kuleuven.be 
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Co RoT space mission (C onvection, Rotation and planetary Tran- 
sits, |Mdlund|er2[200^ see "the CoRoT book") during the Initial 
Run (IRaOl). We present here a detailed study of AU Mon based on 
a long, uninterrupted series of high-precision CoRoT photometry 
obtained in 2008 as well as on high-dispersion echelle spectra se- 
cured at three ground-based observatories in 2007. 

AU Mon (HD 50846, HIP 33237) is an interacting, eclips- 
ing and double-lined spectroscopic binary consisting of a Be star 
and an evolved late-type giant star which in all probability fills its 
Roche lobe and loses matter towards the Be star. To avoid con- 
fusion, we shall hereafter denote the mass-gaining and the mass- 
losing components of AU Mon as the B- and G-stars respectively 
(We will indeed demonstrate that the cool component is a G- 
type star rather than an F-type star which was its classification 
based on photographic spectra). The orbit is circ ular and the or - 
bital period is ll"? 11302 (see O-C diagramJIof iKreinerl I2OO4I) . 
AU Mon is a rare Algol-type system me mber of the W Serpen- 
tis s ubclass, phenomeno logically defined bv lPlavec & KochI l ll978l) 
and lPlavec et"ai] ( Il980h . These binaries are sometimes also called 
"massive or hot Algols". According to current knowledge, they are 
characterized by a semi-detached configuration and several distinct 
components of circumstellar matter: a "cloud" of very hot plasma 
manifesting itself by emission lines in far-UV spectra and probably 
located outside the orbital plane, a disk, possibly optically thick, 
which mimics a false photosphere with a lower r^ff than the central 
star, a gas stream between both components, and a "hot-line region" 
arising fr om the interaction of the gas stream and the disk/cloud 
complex teisikalo & Matsudall2007l) . We mention only som e pub- 
lished studies of W Ser binaries, relevant to our present work. lEliasI 
d 19901) studied centimeter observations of 6 different W Ser bina- 
ries and found very strong evidence for circumstell ar matter around 
these systems. The s ame was claimed for SX Cas l l Andersen et al.l 
Il989l) and RX Cas jAndersen et al.iri989l) on the basis of photo- 
metry and spectroscopy. A system which is simil ar to AU Mon is 
WCru, consisting of a B-star and a G supergiant. IPavlovski et al.l 
( I2OO6I) found this object to have a very extended accretion disk with 
a clumpy rim, and suggested that the dumpiness may account for 
the light curve distortions and asymmetries, as well as for secular 
changes. 



2 PREVIOUS KNOWLEDGE ABOUT AU MON 
2.1 Light changes, orbital period and ephemeris 



The e clipsing nature of AU Mon was disco vered by iHoffmeistei 



1931 



). The first ephemeris was published bv lFlorial (Il937h 



1938 ) derived a slightly longer value of the period. Lorenzil 



Lause 
I977I) 



observed it and concluded that its light curve was undergoing rapid 
cyclic changes. He suggested that this could be due to a fast apsidal 
motion with a period of 243'!2 but also po inted out problems with 
such a n interpretation. Shortly thereafter, ICerruti-Sola & Lorenzil 
( 1 19771) came up with the connect explanation: the brightness of 
the who l e system varied secularly. To understand the phenomenon, 
iLorenzil ( fl980bl) secured 2616 V differential observations between 
1976 and 19 79 at t wo observatories. From a detailed analysis 
of these data, iLoren zi ( 1980a) determined that the brightness of 
AU Mon varied cyclically with a period of 41 1 days and a peak-to- 
peak amplitude of about 0"?2. Taking this variation into account, he 
derived an improved ephemeris for the binary system: 



' |http : / /www ■ as ■ wsp ■ krakow . pi /o-c /index . php3| 



Table 1. Spectros copic orbital elements for AU Mon derived by 
ISahade et al ] il997l) in a spectroscopic analysis. 



Element 



Value 



ag sin i (10 km) 
aa sin((10*km) 
/s:s(kms-l ) 
Kq (km s" 



) 



y (km s"' 
Mb (Mo) 
Mg (Mg) 



6.6 ±0.39 
2.2 ±0.41 
43 ± 2.5 
147 ± 3.0 
0.06 ± 0.02 
204 ± 22 
2.1 ± 1.5 
6.1 ±0.61 
1.8 ±0.53 



TMrni = JD2442801.3752(51) + Il'?1130371(68) x £, (I) 

which has since been used by various other investigators. ILorenzil 
l ll985h analyzed new observations from 1983 and 1984 and 
confirmed the periodic character of the brightness variations of 
AU Mon. A regularly updated ephemeris, based on all compiled 
times of minima, can be found in the database of the O-C 
gramsl iKreineill2004t) '. 

2.2 Spectroscopic investigations 

The first spectroscopic study aim ed at the determination o f the or- 
bital elements was carried out bv lSahade & Cescol ( Il945l) . Helped 
in the spectral classification by Dr. Morgan, they concluded from 
the spectra, at maximum light and during the eclipse, that the pri- 
mary is a B5 main-sequence star and the secondary has a spec- 
tral class near FO. In a later studv. lSahade & Ferreil ( 1198^ investi- 
gated not only ground-based but also the far-UV spectra of the bi- 
nary obtained with the International Ultraviolet Explorer (hereafter 
lUE). Faint line s of the second ary were detected in the Na I 5889 & 
5895 A double t JPoDpeil 1962 ) and later on measured quantitatively 
l |PopDeilll98^ . |PoDpeil (l I962n also noted the presence of a double 
Ha emission. Th e first study, partly based on electronic spectra, 
was published bv lSahade et al.l ir997() . A previous determinations 
of the orbital elements of AU Mon is listed in Table [T] 

The far -UV spectra of AU Mon were obtained by l UE an d 
studied by 



Polidan & PetersI 
1 I 1 



Peters & PolidaiJ (Il982l) : 
Peters (1988); Egikyan (1989); Peters ( 199L 'l994a''bl l 19961) and 
Peters & Polidan (1998.) . Peters & Polidan (1982) confirmed that 



the gas stream between the components is seen in the UV resonance 
lines of N V and C IV in the orbital phase range 0.85 - 0.9 3, where 
phase zero corresponds to the time of primary eclipse. lEgikvanl 
( Il989l) found the electron density of the gas in the envelope to be 
~ 2 - 10" c m~^ , and the electron temperature to be T^ ~ 20 000 K. 
IPetersI flgsj) found a correlation between the instantaneous mass 
transfer rate and the line and continuum spectrum along the 411- 
d cycle. At the maximum of the mass transfer rate, the object 
was fainter and the accretion disk denser while, at the brightness 
maximum, the high-temperature plasma was more prominent in 
the UV spectral lines. Analyzing the spectral energy distrib ution 
(SED) from the low-dispersion lUE spectra, IPetersI ( Il99ll) con- 
cluded that the optical brightness variations of the system are due 
to a 1200 K variation in the photospheric effective temperature of 
the Be primary. She also suggested the presence of a 10 000 K con- 
tinuum source from an op tically thick a ccretion disk around the 
primary to model the SED. IPeter^ d 1994^ lbl) increased the estimate 
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of the temperature of that source to 12 000 K and suggested that the 
cyclic variations in the mass tr ansfer rate could be due to pulsation 
of the mass-losing secondary. IPeters & PolidanI ( Il998h obtained a 
single high-dispersion ORFEUS-SPAS II far-UV spectrum (920 - 
1210 A) and concluded that it is typical of a B3 V star (the SED 
including the lUE data can well be fitted with T^g = 17 000 K and 
log g = 4.0 [cgs]) with V s i n = 1 20 km s"' . 

Richards & Albrightl i l 19991) made an in-depth study of the 



types of accretion structures found in Algol systems. AU Mon was 
part of their Ha spectroscopic study. They have put AU Mon in 
the group with a widely separated double-peaked disk-like struc- 
ture. They have permanent, but variab le, accretion disks s imilar to 
those found in cataclysmic variables. iMiller et"al] | |2005|) studied 
the properties of the accretion structures in AU Mon through a mul- 
tiwavelength spectral study. Double-peaked emission was detected 
in the observed Ha line confirming the presence of an enduring 
accretion disk. The stren gth of the emission varie s with epoch. 

In a careful study, iGlazunova et alj ^20081) derived new ro- 
tational velocities for the components of 23 close detached and 
semi-detached binaries. For AU Mon , they showed how the pres- 
ence of circumstellar matter can falsify the determination of the 
projected rotational velocities. They attempted to avoid the prob- 
lem and obtained i'b sin ; = 124 ±4 kms"', and vq sin ; = 42.1 
±2.1 kms '. They also derived the asynchronicity parameter of the 
B-star, Fb = 5.2. 



2.3 Published light curve solutions 

To date, a few detailed photometric studies of AU Mon have been 
published. Several attempts, all of which are based on Lorenzi's ob- 
servations prewhiten ed for the 41 1-d variation, were made to derive 
light curve solutions. lLorenzil ( ll982h himself made the first attempt, 
creating symmetric normal points from his light curve and using 
the Russ ell-Merrill method. In dependently, and a few months be- 
fore him. lGiuricin et alj ( Il982h solved the light curve using Wood's 
model. Finally, another solution, this time based on the Wilson - 
Devinney (WD) method ( Wilson & Devinnev" 197 lll\yilsonll 19941) , 



was published by Vivekananda Rao & Sarma dl998l) . These au- 
thors concluded that no third light is present in the system. The 
main results of these studies are summarized in Table|2] We note a 
good agreement between the first and the third solution, which are 
based on two independent computer programmes. 



Table 2. Published light curve solutions based on lLorenzi 's [r980bl obser- 
vations. The symbol q = M^jM^ stands for the mass ratio while the sym- 
bols a, b, and c denote the relative dimensions of the triaxial ellipsoids, 
and the back, side, and pole radii for the W D model. The followin g codes 
are used to id entify the authors: GMM82... [ Giuricin et al, h982l) : L82... 
iLorenzil h982l) : VRS98... lw"vekananda Rao & Sarmal il998l) . 



Element 


GGM82 


L82 


VRS98 


in 


78.4±().5 


80 assumed 


78.74±0.06 




0.115±0.O13 


0.18 


0.131±0.001 






0.18 




Clp 


0.267±0.056 




0.2741±0.0034 


bv 


0.242±0.038 




0.2417±0.003 


CF 


0.232±0.031 




0.2324±0.0031 


7-eff (B eq.) (K) 


15 000 assumed 


15500 assumed 


14500±1000 


7-,ff (B pole) (K) 


15010 






7-eff (Feq.) (K) 


6 600±150 


5300 


6000±40 


7-eff (Fpole) (K) 


16860 






q 


0.2 




0.1985 


Lb 


0.645 


0.93 


0.6590±0.0051 


Lf 


0.355 


0.07 


0.3410 



Table 3. The 6 new times of minima for AU Mon derived from the CoRoT 
data. The O-C values were calculated with our new ephemeris (2) for 
AU Mon. 



Time (BJD-2454000) 


O-C 


136.67003(1) 


-0.003 


147.77625(1) 


-0.010 


158.89521(1) 


-0.004 


170.00828(2) 


-0.004 


181.12661(1) 


0.001 


192.23762(2) 


-0.001 



primary minima from the CoRoT light curve. Table|3]lists the 6 new 
CoRoT minima. The published visual, photographic, photoelectric 
and CCD times of minima, reproduced in Table lAll were obtained 
from the General Search Gateway of the Variable-Star Section of 
the Czech Astronomical Societjlj where also references to original 
observers can be found. 



3 THE NEW PHOTOMETRIC DATA AND THEIR 
ANALYSIS 

The new photometric data set consists of a continuous series of 
139704 individual CoRoT photometric observations spanning an 
interval of 56 days, i.e., 5 orbital periodj^ The spectral domain of 
CoRoT enfolds the range from 370 to 950 nm, the average time 
sampling was 32 s and roughly 10% of the datapoints were deleted 
because they were extreme outliers. An estimate of the noise level 
of the light curve, computed as the average of the periodogram be- 
tween 80 and 120d"' is 60/jmag. In addition to the CoRoT data, 
we critically compiled and homogenized all published photoelec- 
tric observations with known dates of observations as well as all 
available times of minima known to us and derived 6 new times of 

^ The CoRoT data of AU Mon are public and can be accessed at 
Ihttp : / / idoc-corot . ias . u-psud . f r/| 



3.1 Time-series analysis and new ephemeris 

Fig. [T] shows the complete set of CoRoT observations vs. HID. 
We converted CoRoT fluxes (F) to magnitudes (m) using in = 
-2.5 log F -I- Co, with Co a calibration constant. We derived Co 
through a comparison between CoRoT magnitudes and visual mag- 
nitudes from literature for all constant stars in the CoRoT field of 
AU Mon. This gives a value of Co = 23.16 ± 0.05 mag. 

First, we verified that Lorenzi's ephemeris (TJ can reconcile 
both the compiled published photometry and the CoRoT data. Con- 
sequently, we used his ephemeris for the initial analyses. In the final 
modelling of photometry and radial velocities (RV hereafter) with 
the PHOEBE programme, described in Sect.|5] we derived the fol- 
lowing improved orbital ephemeris 

^ Ihttp ://var.astro.cz/gsg| 
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Figure 1. The complete CoRoT light curve of AU Mon. The upper panel 
shows an enlarged segment of the curve with very clear rapid non-orbital 
variations. The noise on the data is smaller than the symbol. 



> 




0.4 0.6 
Phase 

Figure 2. Original V-band light curve phased against the period of 
11.1130374 days. 



HJD 2454136.6734(2) + ll"! 1130374(1) x £. 



(2) 



We collected all the available photoelectric observations of 
AU Mon and put them onto a comparable photometric system. The 
journal of all observations is shown in Table|4l the data itself can be 
found in Table lBll We converted t he uvby data to Johnson UB V us- 
ing the transformation derived bv lHarmanec &Bozicl l l200lh . The 
Hipparcos data were transformed to the V magnitude o f the John- 
son sy stem using the transformation formula derived bv lHarmaned 
( Il998h . Fig.|2]shows the phased V light curve based on all the pub- 
lished observations we could colle ct from the literature, including 
the 2698 observations provided bv lLorenzil ( Il980blll985l) . We can 
observe a systematic shift between two extreme states of the max- 
imum light (at the levels of 8.2 and 8.4 mag respectively) as well 
as of the primary minima (at the levels of 9.0 and 9.2 respectively). 
This effect is clearly due to the long-term periodicity reported by 
lLorenzil ( fT980aill985h . 

Subsequently, using subsets of V data sorted into narrow 
orbital-phase bins, as well as all the V data outside of the phase of 

[ rimary minimum, we carried out a period search using the PDM 
Stellingwerd 19781) and the CLEAN jRoberts etaljl987f) methods. 
Both algorithms yielded the same value of Piong = 417 ± 8 days. 
Fig.|3]shows the V magnitude outside minima plotted versus phase 
of the 417-d period. Our linear ephemeris for the total brightness of 
the system reads as follows 



' max.tot.briijhtness 



: HJD2443105.1(±1.4) -)-416'!9(±8J7) x£. 



(3) 



This also means that the CoRoT light curve was collected at a phase 
of total light minimum (in the phase range from 0.46 to 0.59). 

The V-band light curves were inspected at the phases of light 
minima and maxima. The two curves are only shifted for about 
0"J2 with respect to each other (see Fig.|2]l but have the sa me shape 
and a mpUtude, as already concluded by [Cer ruti-Sola & LorenzH 
( Il977h . This was confirmed by tentative light curve solutions for 
both curves, which led to the same values of stellar radii and or- 
bital inclinati on. This result does not support the idea promoted by 
|Petersl ( ll99l1) , that the observed changes could be related to appar- 
ent changes in the T^ff and radius of the B-star's false photosphere. 
If the 417-d variation were caused by some variable third light, then 
the eclipses should become shallower at the 417-d light maxima. 
Similarly, a precession of the orbit would lead to a changing orbital 
inclination. Since none of all these effects is present, we conclude 



8.2 
I 8.3 
> 8.4 

8.5 


_ 0.05 

£ 0.1 
> 

CD 

0.15 

0.2 
-0.7 



g) -0.6 
E 



1 1 1 


1 1 1 1 










u 


1 1 1 





-0.5 



-0.4 



0.2 



0.4 0.6 O.E 
Phase (416.9 d) 



Figure 3. The V magnitude of AU Mon outside minima plotted vs. phase of 
the 417-d period (top). The B-V and U ~B colour changes along the 417-d 
cycle of AU Mon are shown in the middle and bottom panel (see Sect.[8]for 
details). 



that the long-term variation must be caused by a periodic attenua- 
tion of the total light of the binary by some variable circumbinary 
matter The reason why we do not see emission lines coming from 
this circumbinary envelope might be that the envelope is cool and 
dusty. Probably, this huge envelope will also have a small rotational 
speed. Thus its emission will be blended with the stronger emission 
from the circumstellar matter. 
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Table 4. Journal of available photometry of AU Mon with known times of observations. 



Station Time interval No. of Passbands HD of comparison Source 

(HJD-2400000) obs. / check star 



1 


42790.3^2879.3 


278 


V 


50109/50346 


Lorenzi figSOb) 


2 


43023.6^3843.4 


2143 


V 


50109/50346 


Lorenzi fl980b) 


3 


43877.7^3880.9 


195 


V 


50109/50346 


Lorenzi n980b) 


4 


44284.4^4287.4 


13 


B-V 


all-sky 


Kilkennv et al. (1985) 


4 


44288.3^4289.4 


2 


UBV 


all-sky 


Kilkenny et al. f 1985) 


4 


44677.3^4678.3 


2 


UBV 


all-sky 


Kilkenny et al. (1985) 


5 


45254.9^5322.8 


14 


uvby 


50747/50820 


Manfroid et al. (1991) 


2 


45343.4^5790.4 


82 


V 


50109/50346 


Lorenzi (1985) 


6 


47987.4^9056.1 


81 




all-sky 


Perrvman & ESA (1997) 


7 


54136.0-54192.8 


139704 


COR 


all-sky 


this paper, Corot 



Individual observing stations and photometers distinguished by the running numbers in column "Station": 

I... Torino Observatory 1.04-m reflector, EMI 9502 tube; 2... Torino Observatory 0.60-m reflector, EMI 6256S tube; 3... CerroToloio 0.40-m reflector, EMI 6256 tube; 4... South African Astronomical 
Observatory (SAAO) 1.0-m & 0.5-m telescopes, EMI 6256 & 9659 tubes; 5... European Southern Observatory La Silla telescope; 6... Hippairos Satellite; 7... Corot Satellite. 



4 SPECTROSCOPIC ANALYSIS AND RADIAL 
VELOCITIES 

4.1 Initial analysis of the spectra 

New high-dispersion spectroscopic observations were obtained si- 
mukaneously with the CoRoT data over three consecutive orbital 
periods of AU Mon. They were carried out in the framew ork of the 
CoRoT ground-based follow-up programme ( Uvtterhoe ven et al.l 
l2008l) . The da ta we used here consist o f 16 FEROS echelle spec- 
tra (R~48Q00, [Kaufer et al.|[l997ll 19991) from the 2.2-m ESO/MPl 
telescope at LaSilla, Chile, 1 1 well-exposed echelle spectra secured 
with the SOPHIE spectrograph (R~70000, IPerruchot et"ai]|2008l) 
attached to the 1.93-m telescope of the Haute Prov ence Observa- 
tory, and four spectra taken with FOCES (R~65000. |Pfeiffer et all 
1 19981) at Calar Alto Astronomical Observatory. The journal of ob- 
servations, together with adopted radial velocities, is given in Ta- 
blelD 

All data were subjected to the normal reduction process, 
which consists of de-biasing, background subtraction, flat-fielding 
and wavelength calibration. All the reduced spectra were subse- 
quently shifted to the heliocentric frame. Continuum rectification 
and cosmic-spike re moval were carried out man ually using the 
SPEFO programme iHom et alJll996l : ISkodall 19961) . written by the 
late Dr. Jin Horn and until recently being developed by Mr. J. Kr- 
patfl 

To have some guidance before application of more sophisti- 
cated methods, we first derived the RVs via classical measurements. 
Using the programme SPEFO we carefully rectified all red parts of 
the spectra (between 5500 and 6700 A) and cleaned them from cos- 
mics and flaws. The RVs of both components were then measured 
comparing the direct and flipped line profiles. For the G-star, we 
measured Cal 6102.723 A, Fel 6141.730 A, and Fel 6400.000A 
which are all well-defined, unblended, and relatively strong spec- 
tral lines. The r.m.s. errors of the mean RV of these 3 lines ranged 
from 1 to 4 kms"'. The corresponding orbital RV curve is shown 
by black circles in Fig.|4] The only strong and unblended line which 
seems to be related to the B-star in the studied red wavelength range 
is He I 5875.732 A (open circles in Fig.|4l(. For the Ha line (red dots 
in Fig. I4J, the setting was made on the steep wings of the double 



Table 5. Radial velocities measured from the FEROS (EE), SOPHIE (SO) 
and FOCES (EO) data. The RVs of the G-star were derived via 1-D cross- 
correlation of suitable parts of the red wavelength region while the B- 
star RVs were derived near' 4000 A using a 2-D TODCOR-type cross- 
con'elation. The last three RV columns are the measurements carried out 
in SPEEO - see the text for details. The last column denotes the instrument. 



HJD- S/N 


RVb 


RVg 


RVuel 


RVg 


RVh. I. 


2454000 


(km.s"') 




(kms~') 


(kms-') 


(kms-') 



107 
118 
76 
103 
124 



76.6493 95 
76.6697 112 
78.6665 86 
78.6902 86 

103.7806 117 

104.6490 105 

104.8285 98 

105.6557 

105.8796 

106.6320 

106.8208 

107.7008 

108.7038 129 

108.8141 101 

108.8739 159 

109.6744 93 

110.5600 125 

113.3762 72 

114.3554 

115.3513 

118.3702 

119.4596 

121.3731 

125.7931 

127.6520 103 

128.3489 58 

128.6499 

129.3422 

130.3320 

131.3412 

132.3469 



38 
41 
60 
53 
59 
99 



115 
64 
62 
69 
57 



23. 3± 
26.6± 
49.8± 
59.5± 
-6.553± 
-2.108± 
-7.657± 
■14.990± 
■15.247± 
■14.334± 
■16.287± 
■10.859± 
-2.453± 
3.384± 
0.963± 
16.264± 
36.924± 
26.498± 
19.044± 
-2.481± 
■13.448± 
-5.039± 
26.867± 
■37.327± 
■12.026± 
-9.204± 
■13.868± 
■13.363± 
-2.870± 
8.946± 
27.212± 



4.2 
3.1 
4.4 
6.2 
6.939 
10.153 
8.524 
9.559 
8.629 
7.206 
6.911 
4.584 
0.644 
3.488 
2.589 
4.232 
3.668 
3.610 
5.689 
4.473 
7.440 
5.894 
7.542 
7.089 
9.820 
2.343 
7.290 
6.859 
7.891 
1.713 
6.149 



1.36 
3.30 
2.37 
1.57 
1.78 
2.30 



-76.74± 3.60 

-80.23± 4.06 

-129.55± 3.15 

-128.60± 2.84 

55.94± 1.73 

121.79± 1.04 

132.61d 

168.19d 

172.65d 

171.27d 

164.85d 

125.02d 

42.17± 4.11 

28.12± 3.99 
18.01±11.45 

-57.97± 2.87 

-111.92± 2.67 

-65.29± 1.88 

12.74± 1.48 

94.00± 2.63 

150.47± 2.79 

79.44± 3.00 

-95.42± 3.96 

37.10± 2.34 

161.75± 1.50 

174.38± 1.69 

174.36± 2.42 

155.57± 2.06 

98.79± 1.55 

-9.75± 4.66 

-87.30± 1.37 



41.4 
42.7 
78.2 
91.3 
26.9 
0.6 
14.0 
-7.5 
■20.1 
1.9 
■17.3 
■10.3 
2.0 
3.0 
1.8 
23.1 
35.7 
51.9 
90.7 
20.3 
■26.5 
■13.3 
18.2 
7.2 
2.2 
7.5 
-2.2 
-6.6 
2.5 
5.8 
18.5 



-76.3 
-78.7 
-121.2 
-127.0 
56.4 
124.8 
137.3 
175.2 
174.6 
173.4 
168.3 
126.1 
44.6 
30.2 
26.1 
-55.0 
-105.7 
-59.4 
15.2 
92.5 
150.7 
90.6 
-88.1 
41.0 
164.5 
177.8 
177.3 
158.1 
98.4 
-2.3 
-85.9 



13.0 EO 

16.0 EO 

33.1 EO 
43.8 EO 
-2.7 EE 

14.2 EE 
11.2 EE 

1.7 EE 
3.4 EE 
6.1 EE 

9.1 EE 
-24.3 EE 

1.2 EE 
-5.2 EE 
-6.7 EE 

15.1 EE 

27.2 EE 

74.2 so 
61.4 so 
30.7 SO 
-2.8 SO 

-23.9 SO 

41.3 SO 
31.6 EE 

2.0 EE 
- SO 
-10.0 EE 
-16.8 SO 
-28.0 SO 

0.2 SO 

19.4 SO 



Mr. J. Krpata passed away on February 6, 2009 
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Figure 4. The RV curves of the late-type G-spectrum (black dots), B- 
spectrum (open circles) and the Ha emission wings (red dots) measured 
in SPEFO. One can see that the B-spectrum RV, based on the He I 5876 A, 
and that of the Ha emission wings, define similar RV curves which are both 
shifted in phase with respect to the expected RV curve of the G-component. 



Table 6. Trial circular-orbit solutions for the direct RV measurements in 
SPEFO (top). FOTEL ciicular-orbit solutions for the RV measurements 
via cross-correlation of selected segments of the spectra as described in 
Sect. l4.4l and tabulated in Table[5](bottom). The orbital period was kept fixed 
at ll'? 11 30374. The epoch of the primary minimum is in HJD-2454000, 
rm.s. is the r.m.s. error of 1 observation. 



Element 


G-star 


Ho- emission 


He I 




mean RV 


wings 


absorption 


■^min I 


136.668±0.018 


137.94±0.17 


137.46±0.13 


^(km s-') 


157.64±0.16 


32.7±1.7 


38.6±1.4 


y (km s"') 


22.3±1.2 


19.4±2.4 


23.3±2.4 


r.m.s. (kms"') 


6.19 


12.4 


12.0 


Element 


G-star 


B-star 






mean RV 


mean RV 




^minl 


136.684±0.016 


136.972±0.086 




*:(kms-') 


158.51±0.14 


30.41±0.79 




y (kms"') 


18.8±1.1 


12.2±1.1 




r.m.s. (kms"') 


5.47 


5.98 





emission and the measurements were carried out only when a reli- 
able setting was possible. 

The RV curve of the G-star is well defined and appears sinu- 
soidal. To check on the presence of a possible small eccentricity, 
we first derived trial solutions for an elliptical orbit using the pro- 
grammes SPEL (unpubl i shed, w ritten by the late Dr. Jlrf Horn) and 
FOTEL l lHadrava|[l990ll2004j) . We obtained e = 0.045 ± 0.006, 
£j = 93°5 ± 7°2, and Kg = 1 56.71 ± 0.92 km s"' and the test sug- 
gested by ILucv & SweenevI ( Il97 ih indicated that the eccentricity 
is significant. However, we do believe that the orbit is circular and 
the eccentricity is spurious, caused by a small difference between 
the photocentre of the Roche-lobe filling star and its centre of grav- 
ity and/or by the Rossiter - McLaughlin effect (as detailed in Sect. 
6). It was pointed out bv iHarmane 3l l200lL[2003l) that any distur- 



Table 7. Basic properties of AU Mon estimated from spectroscopy (see the 
text for details). The values of v sin (' were derived using the disentangled 
spectra (see Sect. |4.3k the semi-amplitude Kq = 157 - 159 kms"' was 
derived from our spectra and the orbital inclination was assumed to be in 
the range from 77° to 83° (see Sect.|5). 



B-star 



G-star 





15000±2000 K 


5750±250 K 




3.5±0.3 




M 


6.37+2;18 Mo 


1.17+«|9Mo 


R 


7 is+5.77 Tj 

-2.92 ° 




.^synch. 
' eq. 






V sin (' 


116±2kms-' 


43.8±3.5kms-' 


q = Mq/Mb 


0.17 


+0.03 
-0.03 



bance of the sinusoidal shape of the RV curve which is symmetric 
with respect to the systemic velocity must lead to a formally eccen- 
tric orbit with a longitude of periastron of either 90° or 270° . This 
is exactly what we found. Note also that for these orientations of 
the binary orbit, the photometric minima are separated for exactly 
one half of the orbital period so that even very accurate photometry 

cannot provide an ad ditional test. 

Using FOTEL l lHadravd[l990l , l2004ij) . we therefore derived 
a trial circular-orbit solution for the G star, which is compared 
to formal solutions for the RVs of the Ho- emission wings and 
He I 5876 A absorption in Table |6] The epoch of the primary min- 
imum derived from the G-star RVs agrees with that of (photo- 
metric) ephemeris In contrast to this, the RV curve of the 
He I 5876 A absorption and Ha emission are mutually similar but 
shifted in phase with res pect to the clea n anti-phase orbital mo- 
tion of the G-component. IPopped ( Il989h also suggested such be- 
haviour. All this is the situation reminiscent of another W Ser star 
with a similar orbital period of 1Z!9, namely j8 Lyr, where such a 
behaviour is due to the fact that the absorption and emission lines 
originate in bipola r jets emanating from the region of interaction 
of the gas streams jHarmanec et alj 19961) . This may imply that we 
do not see really photospheric lines of the hot mass-gaining star of 
AU Mon. 



4.2 Mass ratio and secondary mass 

To cope with the above problem, we attempted to obtain some esti- 
mate of the mass ratio independent of the B-star's RV curve. Since 
the contact components of semi-detached b inaries usually r otate 
synchronously with their revolution (see e.g. iHarmanedl 1 988h . we 
made this assumption also for A U Mon and used the procedure 
devised by Andersen et alj ( Il989l) as detailed in the Appendix of 
iHarmaned ( Il990l) to obtain an independent estimate of the mass 
ratio iq), component masses and the radius (Rq) of the Roche-lobe 
filling secondary. Its principle lies in the fact that the relative dimen- 
sions of the Roche lobe depend solely on the binary mass ratio q 
while the absolute radius of the spin-orbit synchronized secondary 
is uniquely given by its equatorial rotational velocity, inclination i 
of its rotational axis (assumed to be identical to the inclination of 
the orbit) and by the rotational (=orbital) period. The binary sep- 
aration is given by the third Kepler law so that there is only one 
mass ratio for which the secondary is just filling the corresponding 
Roche lobe. We assumed a circular orbit, A'q = 158 ± 1 km s"' and 
a range of orbital inclinations from 77 to 83° (see Sect.O. We took 
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the uncertainties of thie parameters Kq, v sin i, and Poib. (used in the 
procedure) into account to provide the error estimates of the result- 
ing values of q, Mq, and Rq. The results are listed in Table|7] The 
mass ratio of AU Mon equals q = 0.1 7 ± 0.03. Our results agree 



with 



yivekananda Rao & Sarmal ( 1 19981 . see Table|2l(. ISahade et al.l 



found a mass ratio of q = 0.29 which is much larger com- 
pared to any other study. This is due to their usage of the value 
of Kb, which is probably too optimistic. As Fig.|7]shows, their B- 
spectrum RVs, from phase 0.0 to 0.5, are much more negative than 
any other recorded RVs. Moreover, our separate circular-orbit so- 
lutions for their B- and G-spectrum RVs give systemic velocities of 
1.6 ± 2.1 kms"', and 5.4 ± 3.3 kms"', differing quite substantially 
from the systemic velocity found by us and in all other previous 
studies. 



4.3 Spectra disentangling and spectral type classification 

To separate the spectra of both binary components in an objective 
way, we used the di sentangling procedure developed bv iHadrava 
( 1 19951 . 1 1991 l2004bl . see references therein - Release 2.12.04 of 
the KOREL programme made available to YF). The disentangling 
was applied to different parts of the spectra. T he orbital period was 
fixed at the value derived bv iLorenzil i 1 980 j . see ephemeris ([TJ of 
the present paper) and a circular orbit was assumed. 

Within the parameter range we investigated, the disentangled 
spectra were sufficiently stable to permit to derive the projected ro- 
tation velocity and the effective temperature of both stars. We there- 
fore disentangled the 5500-5700 A, 6125-6275 A, and the 4000- 
4200 A wavelength intervals. Eleven isolated lines were selected in 
the secondary's spectrum and two in the primary spectrum, and we 
estimated the components' v sin / by measuring the position of th e 
first zero of the line profiles Fourier transform jRover et al.ll2002h . 
The values we found are give n in Table[7]and a re in good agree- 
ment with those measured by iGlazunovaetal] j2008l) . who tried 
to avoid the effects linked to the presence of the companion and 
of the circumstellar matter. Though our result for the primary is 
somewhat smaller, the difference is not significant regarding the 
scatter and the small number of lines available in fast rotating early- 
type stars to carry out such kind of analysis. We then compared the 
component spectra to synthetic ones, computed for different effec- 
tive temperatures by means of the ATLAS LTE model atmospheres 
( ICastelli et al] 1 19971; ICastelli & Kuruc3 12003|) and the SYNSPEC 
( iHubenv&Lan^l 19951) programme. This comparison was done in 
the red wavelength ranges for the secondary, while we mainly fo- 
cused on the bluest range for the primary. Since the studied spec- 
tral ranges do not show any strong log(g)-dependence, the surface 
gravity of the secondary (expected to be evolved) was fixed to log g 
= 2.5. 

Furthermore, we assumed a microturbulence of 2kms"' and 
a solar chemical composition. Since we used the hydrogen and he- 
lium stark-broadened line profiles for the B-star, our conclusions 
are in this case not significantly affected by these assumptions. For 
the G-star, spectra for different iron content (i.e. within 0.2 dex of 
the solar abundance) and within 2kms"' of the adopted microtur- 
bulent velocity were computed and analysed. The differences be- 
tween the models with different iron abundances were found to be 
smaller than those obtained by using different effective tempera- 
tures within the estimated error bar of 250 K. Usually a change 
in microturbulence affects more significantly the spectrum. How- 
ever, the weaker lines which are present in the different wavelength 
ranges we have analyzed, and falling in the linear part of the curve 
of growth, are not affected by a change in microturbulence. These 
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Figure 5. Comparison between synthetic spectra computed for different 
effective temperatures and the observed disentangled spectrum (black dots) 
of the cool G-component (upper panel) and the B-component (lower panel). 
The best fitting temperatures for the G and B-star are 5750 K and 15000 K 
respectively. The absorption lines around 4070 A are coming from an AO 
spectrum (extended photosphere) and are best reproduced with a spectrum 
of 10000 K. 



lines were also used to estimate the effective temperature of the 
G-star and a good agreement was found between synthetic and ob- 
served spectra. 

In Fig.[5] one can see the comparison between the disentan- 
gled and synthetic spectra for a selected wavelength range for the 
two components. The luminosity ratios we estimated in different 
wavelength domains are plotted in Fig. |6] while the astrophysical 
parameters we obtained are listed in Table|7] As a byproduct of the 
disentangling procedure we obtained the RVs for both components 
(see Fig.|7j. For comparison matters, we further assumed a non- 
perturbed stellar evolution and a solar-like metallicity. To determine 
the mass, radius and equatorial rotation velocity at orbit synchro- 
nization o f the B-star we made use of the evolutionary tracks cal- 
culated bv lSchalleretai] ( ll992h . We interpolated the mass and ra- 
dius for the effective temperature and surface gravity of the B-star 
throughout these evolutionary tracks. The results of this method 
can be seen in Table |7] The errors on Mb and are unfortunately 
too large to determine the age of the star. Our result firmly estab- 
lishes that the cooler star is a bright G-giant, and not an F-star as 
was tentatively classified in the past. 

From the lower panel of Fig. [5] we see that the agreement be- 
tween the best synthetic and observed spectra for the B-star is not 
perfect. The B-star spectrum is clearly affected by the circumstellar 
matter that acts like an extended photosphere and produces addi- 
tional absorption lines generally found in AO stars (T^g ~ 10000 K; 
spectral lines around 4075 A are coming from the AO star spec- 
trum). 
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Figure 6. The expected wavelength dependence of the relative luminosity 
of the G-star to the total light of the system. We assumed effective tempera- 
tures of 15000 K and 5750 K for the B primary and the G secondary, respec- 
tively. The full line is based on Kurucz' models. The light contribution at 
6150 A was fixed at 0.2. The luminosity ratio was estimated from detailed 
comparisons of the observed and synthetic spectra at several distinct wave- 
length regions and these values ai'e shown by black dots with errors. The 
black square with an error box represents the luminosity ratio in the V band 
based on the PHOEBE solution. 



4.4 Determination of RVs from whole segments of spectra 

Since the motion of the circumstellar matter seems to be phase- 
locked with the orbital motion of the B-star, and the companion is 
a G-star with many weak but sharp spectral lines, the radial velocity 
determination for the G-star is quite straightforward. To avoid sys- 
tematic effects due to the presence of spectral lines belonging to the 
B-star, we concentrated on the 5500-5850 A wavelength region and 
carried out a 1-D cross-correlation of the observations with a syn- 
thetic spectru m computed for the stellar sp ectral classes derived in 
Sect. 14. 31 (See lDavid & Verschuerenl 19951 for the correlation max- 
imum location.). To estimate the accuracy of our measurements, the 
cross-correlation was performed separately on 7 subparts of this 
spectral range, each sub-region providing one set of radial velocity 
measurements. In Table[5]we provide the mean RVs obtained from 
those 7 regions and their r.m.s. errors. 

Due to the orbiting circumstellar matter and the fewer lines 
present in the spectrum of early-type stars, the RVs of the B pri- 
mary are much more difficult to measure. The systematic study 
of the whole available wavelength range shows a large scatter due 
to: 1. the lack of spectral lines, 2. the presence of a third spec- 
trum related to the matter orbiting the primary B-star (and show- 
ing lines generally formed in AO stars) or/and due to bipolar jets 
or disk/stream interactions (see Sect. 14. It . The best wavelength re- 
gion candidate for the measurement of the primary's radial veloc- 
ities is the one that encompassed the He I 4009 and 4026 A lines 
where the effects of the (cooler) third ('AO') spectrum seem to be 
absent. Note that the G-star is more than 3 magnitudes fainter at 
these wavelengths according to Fig. [6] Still, to reduce any system- 
atic effect due to the faint G-star spec trum, we used a TODCOR- 
like procedure jZucker&Mazehll 19941) using templates that reflect 
the spectral types of the components (see Sect. |4.3t . This operation 
was carried out several times using different templates for the pri- 
mary, with the astrophysical parameters falling within 1-cr of the 
estimated values in order to have an estimate of the accuracy of the 
procedure. The results are given in Table |5] and compared to other 
published measurements in Fig|7] 
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Figure 7. Radial velocities obtained in this study by various methods 
(SPEFO, TODCOR, ID CC, KOREL) compared to published values. The 
top panel shows our new and final RVs for sake of clarity, the bottom panel 
shows the comparison. The black line shows the final PHOEBE model for 
the G-star. The residuals for the G-star are pictured in Fig. ll II 



We also derived FOTEL circular-orbit solutions for these 
RVs. The results are shown in Table [6] One can see that the so- 
lution for the G-star agrees reasonably well with that based on the 
directly measured RVs (see Table |6l(. The solution for the B-star 
RVs again shows a phase shift with respect to the expected instant 
of the primary minimum and a significantly lower systemic veloc- 
ity. The K-value for the B-star agrees within the errors with the 
AT-value estimated from our principal analysis (see Table[6j. Nev- 
ertheless, because of the phase- and systematic velocity shift, we 
will not use the RV curve of the B-star in the determination of the 
binary physical properties. 



5 MODELLING OF THE LIGHT AND RV DATA 

In order to derive physical properties from the combined light and 
radial velocity curves, we used the PHOEBE program me, release 
031dev (with pho ebe-gui-cairo. |Prsa & Zwittenl2005h buil t on the 
2003 WD method l lWilson & Devinnevll97ll : lwilsonll990l) . to per- 
form the linearised least-squares analyses with the differential cor- 
rections approach. We assumed a semi-detached system configu- 
ration (using Mode 5) and no third light nor spots were included. 
From previous considerations, we decided to rely only on the spec- 
troscopic estimation of the mass ratio and the RVs of the G-star 
for this modelling. We also fixed the effective temperature of the 
G-star at the value of 5750 K, since this value was reliably derived 
from the disentangled spectra (see Sect. |4.3l l. 

The orbital period of 11'? 113037 was initially fixed. For the 
bolometric albedos, Ab and Aq, and the gravity darkening coef- 
ficient, Gb, we used their theoretical values corresponding to the 
type of atmosphere a nd to the spectral types of both stars (see e.g. 
IVuckovic et al.l2007l) . For the gravity darkening coefficient Go, we 
had to adopt a value of 0.9, i.e. much higher than the theoretical 
value of 0.32, to obtain a satisfactory fit with the data. This value 
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Table 8. Parameters of AU Mon together with their formal 3<t uncertain- 
ties (standard deviations) resulting from the combined PHOEBE solution for 
photometry and RVs. 



Parameter 


Value 


Uncertainty 


HJDo (2450000. +) 


4136.6734 


0.0002 




11.1130374 


0.0000001 


in 


78.6 


0.6 


1 


0.17" 


0.03 


y (km s" ' ) 


17.8 


0.6 


A(Ro) 


42 


1 


e 


0.01 


0.03 


Teff.B iK) 


14595 


900 




5750° 


- 




7.8 


0.9 




2.12 


0.06 


Fb 


4.8 


0.9 


[Lb/(Lb +Lg)]v 


0.82 


0.09 


[Lb /(Lb + Lg)]coR 


0.84 


0.09 


Gb 


1.0* 




Gg 


0.9* 




Ab 


1.0* 




Ag 


0.5* 




xb{V) 


0.458* 


0.001 


XB (CoR) 


0.304* 


0.001 


xg{V) 


0.465* 


0.001 


XG {CoR) 


0.542* 


0.001 


ysW 


0.236* 


0.001 


VB (CoR) 


0.195* 


0.001 


yciV) 


0.243* 


0.001 


ya (CoR) 


0.243* 


0.001 


Roche Radii [in units of orbital separation] 


rg (pole) 


0.12 


0.02 


rg (point) 


0.14 


0.02 


I'B (side) 


0.14 


0.02 


i-B (back) 


0.14 


0.02 


ro (pole) 


0.22 


0.02 


re (point) 


0.31 


0.02 


re (side) 


0.23 


0.02 


re (back) 


0.26 


0.02 


Absolute dimensions 


log gB 


3.78 


0.09 


log go 


2.50 


0.03 


Mboi.B (mag) 


-3.1 


0.6 


Mboi_G (mag) 


-0.25 


0.11 


Mb (Mo) 


7.0 


0.5 


Mg (Mo) 


1.2 


0.3 


Rb (Ko) 


5.6 


0.8 


Rg (Re) 


10.0 


0.8 



' adopted from spectroscopy 
' assumed 



lies very close to the value expected for a radiative envelope. The 
limb darkening coefficients x a nd y of the B- and G-stars were taken 
from lCastelli & Kurucj j2004l) for the Johnson V curve while new 
limb darkening coefficients specifically computed for the CoRoT 
passband, which are now implemented in PHOEBE(rel. 031dev), 
were used. The CoRoT limb darkening coefficients were derived 



from a specific extension of th e tables of van Hamm3 i ll 9931) for the 
CoRoT transmission function l lFridlund et al.l2006h . and computed 
b y Van Ham me (private communication). The method is described 
in Ivan Hamm e (1993), it is a convolution of the outgoing intensi- 
ties with the CoRoT transmission curve for different angles. The 
errors on x(CoR) and y (CoR) coming from this method are 0.001. 
The heliocentric Julian epoch of the primary minimum, HJDo, the 
effective temperature of the B-star, Teif.B, the inclination, /, the di- 
mensionless potential, fig, the fractional luminosity of the B-star, 
Lb, the systematic velocity, y, and the eccentricity, e, were all set 
as adjustable parameters. For the asynchronicity parameter, Fg, we 
used 4 ± 1 since the projected rotational velocity of the primary 
component appears to be a; 4 times the estimated synchronous ve- 
locity (which amounts to 30 + 15 km s"' if we adopt the B-star's ra- 
dius from the solution presented below). We also checked whether 
leaving Fg as a free parameter improved our model. This was not 
the case. 

For practical reasons mainly, we constructed 10 different sub- 
sets from the CoRoT light curve. Given that the CoRoT light curve 
of AU IVIon was largely oversampled, we have split the continuous 
time-series into 2x5 full orbital cycles selecting either even or 
odd data points in time, thus yielding 2x5 different CoRoT data 
subsets. In this way, no information was lost and we were able to 
obtain formal errors based on 10 independent subsets of the data. 
Thus, the 10 CoRoT data subsets, all the compiled V-magnitude 
photoelectric observations and the ID cross-correlation RVs of the 
G-star from Table|5]were used in an iterative procedure to search for 
a consistent model. Of course the CoRoT data are of much higher 
quality compared to the V-band light curve which is reflected in the 
scatter (uncertainties) on the data. 

First, we simultaneously modelled each individual CoRoT 
data subset with the past V-band light curve. In this step we used 
equal weights both for the CoRoT and the V-band data. This rapidly 
converged towards one possible light curve model obtained from 
averaging over the 10 found solutions. The mean parameters (com- 
puted with standard errors) were then adopted as the starting val- 
ues for the next step in the modelling. Next, we simultaneously 
modelled each individual CoRoT data subset and the past V-band 
light curve together with the cross-correlation RVs of the G-star. 
The weights we assigned in this step to the photometry are now 
inversely proportional to the square of the r.m.s. errors for each 
dataset. This means that the CoRoT data were given a much larger 
weight (0.98) than the V-band data (0.02) throughout these calcula- 
tions. The reduced values for each modelling approximated the 
value = 1.4 for the CoRoT subsets, the fit for the V-band light 
curve had a value of » 2.1. Again, an improved light curve and 
RV model was obtained from averaging over the 10 found solu- 
tions. 

Finally, we improved the orbital period and the epoch of pri- 
mary minimum, HJDo, of AU Mon using PHOEBE and our best-fit 
model: we simultaneously fitted all the data with the orbital pe- 
riod and HJDo set as the only free parameters. This yielded a pe- 
riod of 1 r! 1 1 30374( 1), whic h is alm ost the same as the orbital pe- 
riod determined by iLorenzil l ll980d) and which we adopted as a 
fixed parameter in our final modelling. The final value for HJDq is 
HJDo = 2454136.6734(2). We checked whether inclusion of third 
light in the modelling would improve our solution. This was not 
the case. We also verified that, when leaving the mass ratio q as a 
free parameter, its value was fully consistent with our spectroscopic 
estimate of q. 

The mean parameters (with estimated uncertainties) and the 
resulting physical properties of both components are listed in Ta- 
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Figure 8. The final PHOEBE model plotted in red on top of one part of the 
CoRoT light curve (top) and the V-band light curve (bottom), after removal 
of the long-term period of Piong = 417 ± 8 days, phased against the period 
of 1 1. 1 130374 days. Below each light curve the residuals of the model and 
observations are shown. 



ble[8] The true uncertainties will be larger because they depend 
on unknown systematic uncertainties and parameter dependencies 
which were not treated here. The observed V and CoRoT light 
curves and the synthetic model are shown in Fig. [8] There is a very 
good agreement between the CoRoT observations and the best-fit 
model: the mean residual value in all the subsets is of the order of 
0.02 mag (Fig. |9ll, though somewhat larger near the phase of pri- 
mary minimum. We note however that the expected effective tem- 
perature for a normal main-sequence B-star with a mass of 7 solar 
masses would be higher than 15000 K we obtained, in the order 
of 20000-21000 K. The radius of 5.6 solar radii also appears too 
large, even for the mass of 7 solar masses, it should be of the order 
of 4 solar radii. This may indicate that also from photometry one 
measures the outer radius of the optically thick pseudophotosphere, 
rather than the true radius of the B star. Seen roughly equator-on, 
such a pseudophotosphere would also have a lower effective tem- 
perature. Fig.[7]shows the PHOEBE model for the RVs of the G-star. 
Fig.llOl illustrates the configuration of the binary at three different 
orbital phases. 



CHARATERIZATION OF THE 

ROSSITER-MCLAUGHLIN EFFECT OF THE G-STAR 



As already mentioned in Sect. |4.1| the RV curv e of the G-star 
shows a small Ro ssiter-McLaughlin (RM) effect ( lRossiteilll924l ; 
iMcLaughliijl 19241) at phase 0.5. This red-/blue-shifted distortion at 
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Figure 9. The 0-C residuals from the final solution shown for the complete 
CoRoT light curve. 
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Figure 10. Representation of the B and G-star of AU Mon pictured in three 
different phases. The cross in the middle panel denotes the center of mass. 
The system does not show a total eclipse. 



the secondary eclipse is due to selective blocking of the light of the 
rotating star during an eclipse. When the primary star covers the 
blueshifted (redshifted) half of the stellar disk, the integrated light 
of the secondary appears redshifted (blueshifted). Because of this 
selective blocking of the stellar surface during the eclipse, a skewed 
line profile is created. This change in line profile shape results in a 
shift in RV, which in turn results in the redshift-blueshift distortion 
seen during the eclipse. The effect mainly depends on the projected 
rotation velocity of the star, the ratio of stellar radii, the orbital in- 
clination, and the limb darkening. To analyse this effect we have 
subtracted the orbital solution (solid curve in Fig.|7j from the RV 
measurements of the G-star. The orbit-subtracted RV residuals are 
shown i n Fig.fTTI We use d the analytical description of this effect 
given in iGimenej ( |2006|) to simulate the RM effect. The ratio of 
the stellar radii rs/rc the inclination and the radius of the G-star 
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Figure 11. The orbit subtracted RV residuals (dots) showing the RM effect 
at phase 0.5. The solid line is the simulation of the RM effect with the 
parameters given in the text, the dashed line shows the RM effect after an 
increase of the radius of the G-star. 



relative to the size of the orbit re were taken from our final orbital 
solution. The equatorial rotational velocity of the star was set to 
44 ± 4kms"'. The rotational axis of the G-star is assumed to be 
perpendicular to the orbital plane. 

The result of our best fit is seen in Fig.[IT] With our current 
input parameters we cannot fully explain the RM effect. A possi- 
ble explanation is that we are primarily not dealing with rotational 
effects but with the fact that the photocentre and gravity centre of 
the Roche-lobe filli ng star (G-star) are not identical and coinciding. 
IWilson et al. pointed out that this could also cause such dis- 

tortions in the RVs during an eclipse. Another explanation comes 
from the fact that we have circumstellar matter together with mat- 
ter circulating from the G to the B -component which has not been 
taken into account in the calculations. In order to check this, we 
have set the relative radius of the G-star re as a free parameter. The 
result is that we could fit the RM effect if we increase to 0.44A 
(see Fig.ll It. 



7 VARIATIONS ON OTHER TIME SCALES 

As Fig.|9]clearly shows, a systematic pattern of rapid light changes 
which is most pronounced near the phases of primary eclipse was 
detected by CoRoT at every orbital cycle. This is also the rea- 
son why the O-C's calculated from our final ephemeris for these 
minima are significantly larger than their accuracy (see Table[3l(. 
The fast changes near the primary eclipse are reminiscent of what 
IPavlovski etaU | |2006|) found for W Cru (most evident in the U- 
filter). In accordance with them, we propose that the variations are 
due to a non-uniform brightness distribution, probably seated in the 
accretion disk. We also note that the Ha profiles from very similar 
orbital phases but from two different orbital cycles differ from each 
other - see Fig.[T2l For instance, the additional red-shifted absorp- 
tion feature is more pronounced in the first cycle than in the last 



As already mentioned in Sect. |2:2l|Petersl l ll994al lhl) suggested 
that the cyclic variations in the mass transfer rate could be due to 
pulsations of the mass-losing secondary. In order to investigate this 
we performed a frequency analysis on the residuals of the light 
curve after subtracting our final PHOEBE model shown in Fig.[8] 
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Figure 12. A sequence of the Ho- line profiles showing a complicated ve- 
locity field and cycle-to-cycle variations. Data from 3 consecutive orbital 
cycles can be seen. In the left corner of each plot we list the time of mea- 
surement (HJD-2454100), in the right corner we give the orbital phase cal- 
culated from ephemeris j2). 



The periodogram of the residual data set is completely con- 
taminated by the binary orbital frequency and its harmonics. As can 
be seen in Figs[8]and[9l the residuals of the CoRoT data near the 
phase of primary minimum show much larger variations than the 
residuals outside eclipse. This periodicity is due to the fact that the 
primary minima are not completely well fitted, thus introducing a 
strong signal with the orbital frequency of 1 I.I 1303 days. In order 
to avoid this we did not consider the residuals during the phase of 
primary eclipse and recomputed the Fourier Transform (FT). The 
result is shown in Fig.[T3] The orbital frequency and its harmon- 
ics are still present (from to 3 d"') in the top periodogram. The 
frequency peaks in the range of 12 to 17 d"' are due to the CoRoT 
satellite orbital frequency (/oit.coR = 13.97 d~') and are thus in- 
strumental. We sea rched for frequencies using the PDM-method 
l lStellingwerllll97§) and prewhitened the data for the orbital fre- 
quency and 8 other low frequ ency signals (< 0.2 d~') usin g a spline 
fit to the bin means (see e.g. IWaelkens & Rufeneilll983h . This re- 
sulted in a residual light curve without the binary and satellite sig- 
nature. 

The residuals from the PDM-prewhitening procedure were 
analysed using the Short Time Fourier Transform (STFT) with a 
Hamming window l lHarrislll978l) . The result can be seen in Fig.ll4l 
The STFT was calculated using a window width of 10 days, and 
evaluated in 100 equidistant time points. The peaks at 10.4 (/i) 
and 8.3 d"' (/2) are clearly visible as nearly continuous frequency 
bands. A short drop in amplitude is noticeable in the 8.3 d"' fre- 
quency band around day t = 38 d. The power excess at lower fre- 
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Table 9. Frequencies found in the residuals together with their S/N 
ratio (we refer to the text for explanation). The error estimate 
iMontgomerv & O'Donoghud Il999h for the frequencies is ±0.0002 d"'. 
The eiTor on the amplitude is 0.01 mmag. 
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ID 


Frequency 


Amplitude S/N 




[d-'] l^iiHz] 


[mmag] 


/l 


10.4081 120.46 


0.34 14.8 


h 


8.3033 96.10 


0.25 9.2 



quency cannot be attributed to one or several stable frequencies. In- 
stead, the time-dependent behaviour of low frequencies with short- 
living amplitudes are visible, resulting in a smeared out region of 
power excess in a full Scargle periodogram (see middle panel of 
Fig.ll3t. It is difficult to draw firm conclusions about the origin 
of these features, due to the added uncertainties from the PDM- 
prewhitening method in this frequency band. However the two 
frequencies in Table|9] are unaffected by our PDM-prewhitening 
method, and are definitely present in the Scargle periodogram. 
Their characteristics are listed in Table|9l Since these frequencies 
are not detectable in the spectra, we cannot directly assign a phys- 
ical origin to them. The frequencies /i and fi are likely to origi- 
nate from the B-star because it produces ~ 82% of the light and it 
concerns isolated frequency peaks. Altho ugh, theory does n ot pre- 
dict such pulsations for B-stars (see e.g. iMiglio et al.ll2007h . such 
frequencie s have also been detected in CoRoT data of the Be star 
HD 49330 ( iHuat et alj2"009l) . This result puts AU Mon in the neigh- 
bourhood of the sample of semi-detached Algol-type eclipsing bi- 
naries with an oscillating mass accreting component (oEA stars , 
with RZ Cas as the best studied object iMkrtichian et alj|2007h . 
AU Mon can not be fully classified as an oEA star, because oEA 
stars are at the end of the mass transfer regime which is still not the 
case for AU Mon. 

On the other hand, theory predict s that the G-star should 
produce solar like oscillations (see e.g. lHoudek&Goughll2002l : 
ISamadi et alj2007l) . Since we know the mass, radius and T^s of the 
star, we predict y^ax = MG/Rl-(T^ff.G/5m)-°-^ -3050 ~ 34a/Hz(= 
2.94 d~') to be the frequen cy at which to expect such oscillations 
jKieldsen & Bedding|l99i) . This value is in agreement with the lo- 
cation of the power excess in the Scargle periodogram of AU Mon 
(middle panel of Fig.fTSt. The amplitude of the power excess is not 
in disagreement with the amplitudes of the d iscovered solar-like os- 
cillations in the CoRoT data of red giants JPe Ridder et al. 1 1200^ 
taking into account the flux contribution of the G-star. Due to a lack 
of mode identification and uncertainty on the origin of the oscilla- 
tions, we are unable to exploit the detected frequencies seismically 
without additional spectroscopic information about them. 



8 UNDERSTANDING THE LONG-TERM LIGHT 
VARIATIONS 

Other cases than AU Mon exist in which the discovery of long- 
term periodic or at least cyclic brightness variations (reminiscent 
of AU Mon ) have been claimed, namely those cases where sys- 
tematic photometry in a calibrated photometric system has been 
carried out or where interacting binaries with Balmer emission 
lines and hot mass-gaining components are di scussed. Pro bably the 
first reports of such a beha viour are RX Cas l lKalvlll979l) . TV Cas 
= HP 1486 l lWalteilll979h and V505 Sgr =BD-14°5578 jWalteJ 
Il98lh . However, the latter two systems exhibit only very faint 
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Figure 13. The periodogram of the residuals out of primary ecUpse together 
with the spectral window (top+inset). The middle panel shows the peri- 
odogram after prewhitening for 9 low-frequency signals. The bottom panel 
shows the final residual periodogram after prewhitening all significant fre- 
quency peaks. 
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Figure 14. The Short Time Fourier Transform of the residuals of the CoRoT 
light curve after prewhitening for the orbital frequency. The units of the 
colour scale are in mmag. See text for details. 



single-peaked Ha emission and their long-term variations are prob- 
ably rela ted to the presence of distant third componen ts in these 
systems jRichards & Albrighlll 19991 : 1 Vesper et ai]|200lh . 

In TablellOl we provide basic information on galactic 
emission-line systems with cyclic long-term variability known to 
us. During the past five years, similar systems were also dis- 
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Table 10. Galactic hot emission-line binaries with cyclic long-term bright- 
ness changes. The tabulated mass ratio is always the ratio of the mass of the 
lighter (usually Roche-lobe filling) component to that of the B-star. 



Name BD Pa^b. P[ons Mass ratio Ref 



RXCas +67°244 32?312 516'?06 0.3()±0.05 1,2,3,4 

AUMon -01°1449 ll-JlH 41T?0 0.17±0.03 5 

CXDra +52°2280 6'?696 130''-180'' 0.23 6,7,8 

ySLyr +33°3223 12^94 282'?37 0.225 8,9,10 

V360Lac -i-41°4623 10?085 322'?24 0.163 11,12 



References to original studies identified by tile running numbers in column "Ref.": 
l... lKalviri97^ : 2... Kriz et Jl Ugsdl) : 3... lAndersen et alitTosl ; 4...!Pustvlnik et arf2()07ll: 5... 
this study; 6...[Koubsky el al, <I98C| ): 7... lHorn el alll I992t ; 8... [Richards cl al. i 2000); 9... 
I Harmanecet 111119961 ); 10... lHarmanell2002l) ll... lHill et allll997l ); 12... lLinnell eraill200d l 



covered and studied rather inte n sively in th e Magellanic Clouds 
jMennickent et al .I2003ll2005bl 'a'. 2006, '20#. In spite of all effort, 
the true cause of the long-term cyclic changes remains unexpl ained, 
althou gh several different interpretations were put forward: iKalvl 
( Il979h interpreted the 516- d periodicity of RX Ca s as pulsation of 
the Roche-lobe filling star. iHarmanec et alj ( 1 19961) considered the 
possibility that the 282-d cycle of j8 Lyr is a beat period between the 
orbital period and rapid changes with a cycle of 4^1 which they de- 
tected in spectroscopy. They also pointed out some similarity with 
the 164-d cycle of V13 43 Aql = SS 433, a massive 13^.08 X-ray 
binary with bipolar iets. lWilson & van Hammd [l99^ also studied 
P Lyr. They concluded that neither apsidal advance nor precession 
can account for the 282-d light variation. They were unable to ex- 
clude pulsations of the disk as the cause. As already mentioned, 
IPetersI ( Il99lh concluded that the optical brightness variations of 
AU Mon are due to a 1200 K variation in the photospheric effec- 
tive temperature of the B prim ary. 

iMennickent et alj j2003h reported the discovery of eclipsing 
binaries with long-term periodic bright ness changes from OGLE 
photometry of the Magellanic Clouds. IMennickent et al.l ^20031) 
pointed out that the light maxima of the long cycles are al- 
ways acconi£anied_by a mild reddening of the objects in ques- 
tion. IMennickent et al.l ( l2005bl) concluded that the long brightness 
changes a re often cyclic rather tha n strictly periodic ones. Sub- 
sequently, IMennickent et al. I (|2008l) concluded that the long-term 
variations must be due to variations in circumbinary matter as we 
found for AU Mon. They suggested that the system experiences su- 
percycles of mass outflow which lead to replenishment of the cir- 
cumbinary envelope. 

AU Mon may become a key object to study the true nature of 
the long-term brightness changes since it is bright and the ampli- 
tude of the 417-d period is large. In spite of our effort to collect 
and homogenize existing photometry, available material on colour 
variations is, unfortunately, very scarce. In Fig.|3]we compare the 
V magnitude and B—V and U - B colour changes along the 417-d 
cycle for all observations we were able to transform into compara- 
ble UBV magnitudes. Unfortunately, there is no colour information 
near phases of the brightness maximum. A plot of both colour in- 
dices vs. time shows that there could also be colour variations on 
longer time scales, than those related to the 417-d cycle. If this 
is a more general pattern, this could perhaps also explain some- 
what contradictory reports on the colour behaviour of various ob- 
jects studied by Mennickent and his collaborators. Clearly, future 



systematic calibrated multicolour photometry of AU Mon over the 
whole 417-d cycle is very desirable. 

We tentatively suggest that the long-term brightness changes 
of AU Mon must be associated with some circumbinary matter and 
it is plausible to assume that the bulk of such material is associated 
with putative bipolar jets for which we find indirect support from 
the shi ft of the RV curves of some spectral lines associated with the 
B-star. iBisikalo et al] ( 120001) showed via hydro calculations for /3 
Lyr that when the encircling stream hits the denser primary stream 
from the mass-losing star, it gets ben ded and goes out of the orbital 
plane resulting in jet-like structures. IPetersI ( |2007|) reported proba- 
ble detections of bipolar jets of very hot plasma, perpendicular to 
the orbital plane, for three other hot interactin g binar ies: TT Hya, 
V356 Sgr and RY Per. However, iMiller et alj ( |2007|) studied cir- 
cumstellar matter of TT Hya and argued against jets. The variations 
themselves could have two possible causes: either cyclic ch anges in 
the ma ss outflow from the binary as suggested by Mennicke nt et al.l 
( l2008l) or precession of the binary orbit which would change the at- 
tenuation of the binary due to changing projection effects of jets. 
This latter idea seems improbable, however, in the light of our 
finding that the light curve from the maxima and minima of the 
417-d period lead to the same binary elements including the orbital 
inclination. It is clear that continuing systematic observations of 
AU Mon, including spectro-interferometry and polarimetry, could 
help to understand the nature of the remarkable changes of these 
interesting objects. 



9 SUMMARY 

Our analyses of very accurate CoRoT space photometry, past John- 
son V photoelectric photometry and high-resolution echelle spec- 
tra led to the determination of improved fundamental stellar prop- 
erties of both components of the massive and interacting system 
AU Mon. We derived new and accurate ephemerides for both the 
orbital motion (with a period of I I'll 130374) and the long-term, 
overall brightness variation (with a period of 4I6'!9). It is shown 
that this long-term variation must be due to attenuation of the total 
light by some variable circumbinary material. We derived the bi- 
nary mass ratio Mq/Mb = 0.17±0.03. Using this value of the mass 
ratio as well as the radial velocities of the G-star, we obtained a con- 
sistent and coherent light curve solution and a new estimate of the 
stellar masses, radii, luminosities and the effective temperatures. 

We must point out that our final model does not include the gas 
stream and accreting matter on the B-star. It would be interesting to 
consider such complications in any future modelling of the binary. 
We also report the discovery of rapid and periodic light changes vis- 
ible in the high-quality residual CoRoT light curves. The rapid light 
changes visible in the residuals near primary minima repeat at every 
orbital period. They are probably due to a non-uniform brightness 
distribution, seated in the accretion disk. Outside the primary min- 
ima of the CoRoT light curve we detect two frequencies, they are 
in the expected frequency domain of B -stars. Complementary in- 
terferometric and polarimetric observations will be needed to even 
better understand the geometry and the nature of the circumbinary 
matter in AU Mon. 
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APPENDIX A: TIMES OF MINIMA 

APPENDIX B: AVAILABLE PHOTOELECTRIC DATA 

SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online ver- 
sion of this article: 

Table lBll Available photoelectric observations of AU Mon. 
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Table Al. Visual, photographic, photoelectric and CCD times of minima. 
The times were adopted form the General Search Gateway of the Variable- 
Star Section of the Czech Astronomical Society (http://var.astro.cz/gsg). 
The 0-C values were calculated with our new ephemeris (2) for AU Mon. 
The third column lists the type of observation, vis: visual; pg: photographic; 
pe; photo-electric. 



liinc 




lype 


1 ime 




lype 


(_rij U-Z^UUUUU ) 












26743.160 


0.124 


vis 


37378.364 


0.151 


pg 


26743.390 


0.354 


pg 


37400.351 


-0.088 


pg 


26765.320 


0.058 


pg 


37578.633 


0.385 


pg 


26976.400 


-0.010 


vis 


37667.505 


0.353 


pg 


26987.530 


0.007 


vis 


38322.601 


-0.220 


pg 


27076.430 


0.003 


vis 


38378.568 


0.182 


pg 


27120.890 


0.011 


vis 


38411.508 


-0.217 


pg 


27487.583 


-0.027 


vis 


38467.353 


0.062 


pg 


27843.110 


-0.117 


vis 


38856.363 


0.116 


pg 


27887.631 


-0.048 
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39200.388 


-0.363 


pg 
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0.099 


pg 


39945.411 
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pg 


28209.912 
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-0.141 


pg 
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-0.143 
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37367.333 


0.233 
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51513.5290 


-0.468 


ccd 



Table Bl. Excerpt of the available photoelectric observations of AU Mon. 
The full Table is available as Supporting Information. 



Data from Lorenzi ^98^ 



Time 
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(HJD-2400000) 


(mag) 


45343.4114 


8.380 


45346.3339 


9.137 


45346.4262 
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45351.4059 


8.415 


45351.4202 
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45353.3960 
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45354.4253 


8.387 


J J J J . joo / 


8 385 


45356.3518 


8.448 


45357.2635 


9.118 


45357.2823 


9.151 


45357.3028 


9.165 


45357.3205 


9.181 


45357.3349 


9.180 


45357.3471 


9.187 


45357.3600 


9.188 


45357.3739 


9.185 


45357.3927 


9.192 


45357.4045 


9.180 


45357.4161 


9.169 


45357.4287 


9.165 


45357.4403 


9.155 


45357.4565 


9.147 



